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ABSTRACT 



The Normalized Radar Cross-Section (NRCS), the fundamental measurement made 
by radar scatterometers, was obtained as part of the Water-Air Vertical Exchanges 1987 
(WAVES87) experiment. The experiment was designed to evaluate the effects of 
environmental parameters on the NRCS and was performed from a research tower located 
in Lake Ontario, on which tw'o microwave scatterometers operating at 14.0 and 5.0GHz 
were installed for six weeks in the autumn of 1987. The novel aspect of this experiment 
was that the 14.0GHz radar automatically rotated through 300° in azimuth angle at six 
different incidence angles to the water surface, accompanied by simultaneous 
measurements of wind stress and high resolution directional wave spectra. Therefore, the 
incidence and azimuthal angle behavior of the NRCS was examined as a function of wind 
speed, friction velocity, wind direction, wave direction and atmospheric stability. 

The dependence of the NRCS on wind speed for various incidence angles is similar to 
previous results. However, the slope exponents of the NRCS vs. 19.5m wind speed 
curves at intermediate incidence angles are higher than the corresponding open ocean 
measurements. Scaling the lake neutral wind speed data by the ratio of lake to ocean drag 
coefficients reduces the slopes of the curves and suggests the drag coefficient has a sea 
state dependence. The correlation between NRCS and neutral wind speed at lm is higher 
(0.91) than between the NRCS and friction velocity (0.73 at 40°). The minima in the 
sinusoidal modulation of the NRCS as a function of relative wind angle (the angle between 
the wind and antenna directions) are often shifted (by as much as 45°) such that the minima 
do not always occur at cross-wind angles. Instead, the angular distance between the 
NRCS minima in the case of a wind-wave sea appears to approximate the directional spread 
of the waves about the upwind direction, generally rather less than 180°. The degree of 
sinusoidal modulation of the NRCS with relative wind angle is highly correlated with 
significant slope and inverse wave age at 20° incidence angle (0.90) and moderately 
correlated at 40° (0.75): i.e., increased azimuthal modulation at 20° is associated with a 
steeper wave field. The dependence of the NRCS on atmospheric stability shows the 
NRCS to decrease by about 5dB between air- water temperature differences of about -16 to 
+ 10°C. This stability 7 effect is removed by parameterization of the NRCS in terms of either 
the friction velocity or neutral wind speed at lm, with the neutral wind speed providing the 
best normalization of the data. The results show that radar scatterometers are an especially 
sensitive means by which to study the air-sea interface: the magnitudes of the 5.0GHz and 
1 4.0GHz NRCS respond nearly instantaneously to changes in the near-surface neutral 
wind speed, but the directionality of the (Ku-band) NRCS is the result of complicated 
interrelationships among the influencing environmental variables. 
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I. INTRODUCTION 



A. REMOTE SENSING OF THE OCEAN 

In the three decades since the launch of the first meteorological satellite, 
the enormous potential of remote sensing of the Earth's environment has 
been proved. Products of visible and infrared imagery are routinely used to 
monitor the evolution of weather systems, oceanic fronts and eddies. Passive 
infrared and microwave radiometers provide vertical temperature soundings 
of the atmosphere and information on its moisture content. Active 
microwave radars serve the dual purpose of measuring sea surface elevation 
and waves, and inferring sea surface winds over areal expanses, both with 
accuracies never before realized. To capitalize on these now proven 
techniques, space agencies have initiated ambitious programs to build and 
launch environment observing satellites in the early 1990's. 

At the heart of these programs are intriguing and unresolved 
questions, answers to which will provide the basis for next generation 
instrumentation and delimit the applicability of remotely sensed data. 
Questions raised are 1) theoretical: How are the spectral emissions/reflections 
from the earth related to geophysical variables of interest?, 2) technical: How 
can advances in computer and engineering sciences be incorporated into 
remote sensing instruments?, and 3) operational: How are the huge volumes 
of data to be received processed, distributed and analyzed in near real-time? 
The issue addressed in this study is of the first general type. In particular, the 
purpose of this study is to evaluate environmental effects on the functional 
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form which relates radar echoes from the sea surface to the sea surface wind 
vector. 

B. THE RESEARCH PROBLEM 
1. Statement of the Problem 

Radar backscatter at intermediate incidence angles (20 to 70 deg) can 
be related by a "model function" to the sea surface wind speed and direction 
due to backscatter which is produced primarily by "Bragg resonant" reflections 
of the incident radiation from successive crests of short, wind generated 
capillary-gravity waves (0.2 to 20 cm wavelength). For the signal to be 
scattered to the radar, these waves must have components traveling along the 
direction of the antenna line-of-sight. Therefore, the backscattered power 
displays a strong angular dependence: maxima occur when the antenna is 
pointed either upwind or downwind, minima occur when the antenna is 
pointed crosswind. Furthermore, the upwind maximum is usually larger 
than the downwind maximum; this difference is enhanced for horizontally 
polarized radiation. The idea, then, behind scatterometer wind vector 
determination, is that a particular area on the ocean surface viewed from 
several different directions; i.e., azimuth angles, yields different wind 
components per azimuth that can be used to determine a single wind vector. 

The above information has been incorporated into previous model 
functions by writing the normalized radar cross-section (NRCS) of the sea as a 
function of radar incidence angle, radar polarization, and azimuth angle 
relative to the true wind vector (Jones, et al., 1977; Moore and Fung, 1979; 
Schroeder, et al., 1982). However, Daley, et al., (1984) identify the deficiencies 
in this simple type of formulation and list the additional dependencies which 
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may need to be included in the scatterometer model function. Based on 

recent research, these dependencies are: 

(1) cross-section dependence on wave slope and wave direction relative to 
the wind vector (Keller, et al., 1985, 1989; Li, et al., 1989), 

(2) cross-section dependence on air-sea temperature difference, and thus 
on stability of the marine surface layer (e.g., Keller, et al., 1989), 

(3) a lack of self-consistency of the model between vertical and horizontal 
polarizations (Wentz, 1984; Woiceshyn, et al., 1986), 

(4) decreased sensitivity of the model function at low wind speeds (below 7 
m/s), (Woiceshyn, et al., 1986), 

(5) cross-section dependence on surface contaminants and films 
(Huhnerfuss, et al., 1983), and 

(6) cross-section dependence on viscosity of the water as a function of 
temperature (Woiceshyn, 1986; Donelan and Pierson, 1987). 

For the scatterometer wind velocity to meet the stringent operational 

requirements of the Navy (±2m/s, +2.5°) , an updated model function is 

required. 

2. Research Objective 

The effects of the various environmental parameters listed above are 
represented schematically in Figure 1. An increase in wind speed produces an 
associated increase in the NRCS through amplification of the capillary-gravity 
waves from which the electromagnetic waves are reflected. Observations 
show that the wind speed dependence can be represented as a power law 
(Figure 2a). Referring to Figure 2b, the cross-section varies with wind 
direction as a result of the anisotropy of the capillary-gravity waves with the 
relative maxima occurring when the radar beam is directed perpendicular to 
the wavecrest (upwind/downwind) and the minima occurring when the 
radar beam is parallel to the wavecrests (crosswind). This behavior is rather 
successfully modeled as a truncated Fourier cosine series. 
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Figure 1. Environmental Effects on the Normalized Radar Cross-Section 

If the short waves are tilted towards the radar and modulated in 
amplitude and frequency by the slope of the longer waves on which they ride, 
the NRCS will again be increased (Figure 2c). Observational evidence for the 
long wave influence has been reported by Keller, et al., (1985,1989) who show 
that the cross-sections increase with wave slope only when the atmosphere is 
stably stratified. Measurements during neutral and unstable conditions 
yielded no discernable dependence of the NRCS on long wave slope; it is 
suspected that the influence of stability is a consequence of the coupling 
between the surface waves and wind turbulence (see Figure 2d). 
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Figure 2. Functional Dependence of Normalized Radar Cross-Section on a) 
Wind Speed for Various Angles of Incidence, 0 (Schroeder, et al., 1982); b) 
Wind Direction (Jones, et al. 1977); c) Long Wave Slope (Keller, et al., 1985); 
d) Atmospheric Stratification (Keller, et al., 1989) 
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Unstable atmospheric stratification, as indicated by a negative air-sea 
temperature difference (Figure 2d), results in 1) a slightly higher near-surface 
wind speed, 2) a wind profile which deviates from the neutral stratification, 
logarithmic form, and 3) an increase in the surface drag coefficient. The 
NRCS is seen to decrease in response to these changes in the surface layer as 
the stratification changes from unstable to stable. 

The dissipation and growth rates of the capillary-gravity waves have 
been shown to be proportional to water viscosity (Lamb, 1932) which is an 
inverse function of water temperature; therefore, this mechanism has been 
proposed to explain observations which show an increase in the radar 
backscatter for increasing surface temperature, particularly at low wind speeds 
(Woiceshyn, et al, 1986; Donelan and Pierson, 1987; Kahma and Donelan, 
1988). Keller and Plant (1988) present recent wavetank measurements which 
do not show a pronounced temperature dependence. Additional 
observations are needed to reach a definite conclusion. 

The physical processes which influence radar backscatter all fall into 
the realm of describing, modeling and measuring the microphysics of the 
air/ water surface layers. Among other institutions which study these 
processes, the Canada Centre for Inland Waters has supported air-water 
interaction experiments since 1976 under its WAVES (Water-Air Vertical 
Exchange Studies) program (Tsanis and Donelan, 1987). The primary 
platforms for the CCIW program are laboratory wavetanks and a tower at the 
western end of Lake Ontario, positioned 1.1km offshore from Hamilton, 
Ontario. In the autumn of three consecutive years (October-December, 1985- 
1987), the tower was outfitted with instruments to study wave spectra, wind 
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stress, near-surface wind profiles and water turbulence. The Naval Research 
Laboratory added two radar scatterometers (14Ghz and 5 Ghz) as elements of 
the 1986 and 1987 experiments in order to take advantage of the opportunity 
to obtain collocated radar and environmental measurements. In 1986, both 
radars were mounted on the tower railing at fixed incidence and azimuth 
angles. In 1987, the 14Ghz radar was mounted on an antenna rotator and 
extended on a boom away from the tower. This configuration allowed the 
radar to be rotated in azimuth at different incidence angle settings. 

The objective of this study is to examine the effects of the parameters 
shown in Figure 1 on measured scatterometer returns from an analysis of the 
WAVES data and to determine if these effects are successfully predicted by 
currently available model functions. Explanations are offered for the observed 
and predicted behavior of the radar cross-section. 

3. Limitations of the Study 

It should be noted that the above list of parameters which affect the 
cross-section is not complete (Masu'.o, et al., 1986). In effect, any oceanic 
phenomenon which has a surface expression may alter the cross-section; e.g., 
surface slicks caused by collection of surfactants by internal waves or 
Langmuir circulations (Huhnerfuss, 1983; Hughes and Gower, 1983), surface 
currents (Phillips, 1981), and the surface adjustment due to topographic 
interactions in shallow coastal areas (Valenzuela, 1985). Similarly, 
atmospheric phenomena such as planetary boundary layer rolls and or coastal 
fronts can be equally important in generating spatially varying surface 
signatures (Geernaert, 1990). 
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Due to the location of the study site in Lake Ontario, the atmospheric 
and oceanic variations described in the above paragraph will not be addressed 
in this study. The maximum wind speed observed, 15 m/s, is not as great or 
sustained as would be possible in the open ocean. On Lake Ontario, the 
longest wave period observed is 8s and "wind rows" aligned with the wind, 
believed to be due to Langmuir circulations, are the most common form of 
surface slick. The fetch limits are the geographic bounds of the lake: 1 to 
300km. The wind direction is primarily from the west; i.e., the short fetch 
direction. Finally, the water is fresh versus saline. At the radar frequencies 
used in the experiment, the dielectric constants of fresh and saltwater are very 
nearly the same, therefore, the cross-section will not be affected. However, 
residence time of foam in a breaking wave regime is longer in saline water, 
which could cause the backscatter observed in the ocean to differ from that on 
the lake, especially at higher incidence angles and wind speeds (Gucinski, 
1986). 

C. NAVAL RELEVANCE: REQUIREMENTS FOR THE WIND VELOCITY 

AT THE SEA SURFACE 

Marine wind is the only environmental parameter listed as a 
requirement in all Naval Warfare areas (Brown and McCandless, 1988). In 
terms of criticality, only sea-state is more important although the correct 
specification or prediction of the sea-state has accurate marine winds as a 
necessary first requirement (Janssen, et al., 1989). Both strategic and tactical 
naval operations are affected by marine winds. c 

The surface wind is an important parameter in strategic battle 
management. For example, the general weather picture (e.g., winds, storms, 
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rain) is essential for decision making, resource management and optimum 
track ship routing. The systematic observation of ocean areas for effective sea- 
control and surveillance is hindered by high winds and waves which degrade 
sonar array performance by increasing ambient noise and reduce radar 
effectiveness by increasing sea clutter. 

All operations at the tactical level can be suspended on account of 
extreme wind and sea-state conditions. In the traditional sense, these tactical 
operations or decisions require only a local measure of environmental 
conditions. However, today’s battlegroup size is 0(1000 km) which extends 
the battle horizon to ocean basin scales when the deployment range of 
modern weaponry is considered. Therefore, next-generation tactical decision 
aids will have expanded areas of coverage. This concept is already being 
implemented in the Tactical Environmental Support System (TESS). 
Shipboard users of this system will use direct readout from overhead satellite 
passes and conventional data to produce local data products. These products 
will be supplemented by those produced at central site locations using all 
available satellite data, model outputs and conventional data. 

The pervasive requirement for surface wind coverage at spatial scales of 
tens to thousands of kilometers argues for global satellite measurements of 
surface winds. An Operational Requirement (OR) for this capability, which 
listed the performance goals shown in Table I (Daley, et al., 1984), was issued 
by the Chief of Naval Operations in 1977. This OR has since been rescinded in 
order to consolidate requirements for environmental parameters, but the 
performance goals and requirements are still valid. The 
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TABLE I. PERFORMANCE GOALS FOR SATELLITE MEASUREMENT OF 

SURFACE WINDS 

NOTE: The double entries indicate the most desirable values on the left and the minimum 



acceptable thresholds on the right of the slashes. 



Surface 


Data Receipt Rate 


Horizontal 


Measurement 


Absolute 


Range 


Wind 


Freq Timeliness 


Resolution 


Precision 


Accuracy 




. Speed 


3//12 hrs 0.25/ /3 hrs 


10//25 km 


5//20 % 


2mps//4mps 


l-75mps 

//3-25mps 


Direction 


3//12 hrs 0.25/ /3hrs 


10/ /50km 


157/10° 


±25°//10° 


0-360° 



most recent compilation of OR's shows at least 24 requirements which 
specify surface wind as an important or critical element. 

D. SUMMARY 

In this study, radar backscatter as a function of wind speed and direction, 
long-wave slope, atmospheric stratification and water temperature is 
examined using data acquired from an instrumented tower in Lake Ontario. 
A historical review of scatterometry is presented in the next chapter along 
with background material on wave and air-sea interaction theory which is 
essential for the interpretation of the results. A review of current model 
function theory and supporting observations is given in the final section of 
Chapter II. A description of the experiment and data processing follows in 
Chapters III and IV. The major results of the study are presented in Chapter 
V. These findings are summarized and the major conclusions stated in 
Chapter VI. 
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IL LITERATURE REVIEW 



A. INTRODUCTION 

To understand the physics of scatterometry, one must couple 
electromagnetic, ocean wave and air-sea interaction theories. It is the purpose 
of this chapter to summarize the aspects of each of these areas which are 
pertinent to scatterometry. The intent is to make clear how certain 
environmental parameters can (theoretically) influence the radar 
measurements and how currently proposed model functions incorporate 
these parameters. Each of the theory sections (C,D,E) derive, in outline form, 
the expressions used in the model function section (Section F). A synopsis of 
the history of scatterometry is presented first in order to provide a larger 
context in which to set the model function problem. 

B. HISTORICAL REVIEW OF SCATTEROMETRY 

In his experiment to test Maxwell’s electromagnetic theory, Hertz, in 1886 
demonstrated that radio waves were reflected from various solid objects. In 
1904 Hulsmeyer (Ulaby et al., 1981) applied this result to the tactical problem 
of ship detection and in the war years to follow, ground-based radars were 
developed to detect aircraft and ships. These radars operated at meter or 
decimeter wavelengths, much longer that those of present-day radar 
scatterometers. To preserve the secrecy of the early radar development 
program, arbitrary letter designations, which are still in common usage, were 
assigned to the microwave frequency bands (Table II). 
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TABLE II. MICROWAVE FREQUENCY BAND DESIGNATIONS 



Letter 


Frequency Band GHz 


Central Wavelength Cm 


p 


.225 


.390 


115 


L 


.390 


- 1.55 


20 


S 


1.55 


- 4.20 


10 


C 


4.20 


- 5.75 


5 


X 


5.75 


- 10.90 


3 


Ku 


10.9 


- 18.0 


2 


K 


18. 


- 26.5 


1.35 


Ka 


26.5 


- 36 


1.0 


Q 


36 


- 46 


.73 


V 


46 


- 56 


.58 


W 


56 


- 100 


.38 



By World War II, airborne radars were used to detect other aircraft 
and ships-at-sea, and some in the 3 and 10 GHz bands were producing images 
of the ground. As the instruments were refined, it became clear that radar 
reflections from ocean waves, which obscured the echoes from the actual 
targets (sea clutter), would have to be accounted for to improve data 
interpretation. Thus, the radars were directed at the sea surface itself with the 
purposes of improving radar performance, developing techniques to study 
the sea surface, and testing new theories which attempted to explain scatter 
from the sea (Stewart, 1985). 

In one of the first controlled radar measurement experiments, 
Goldstein (1946) introduced the standard unit of measurement for radar 
backscatter from the sea, a 0 , the (area) normalized radar cross-section (NRCS). 
Goldstein's results suggested that o° varied with radar frequency, polarization 
and incidence angle (Huebner, et al., 1975). From 1959 to 1975, the Naval 
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Research Laboratory (NRL) continued the efforts to establish the character of 
the scattering cross-section via an extensive measurements program which 
used airborne multifrequency radars to examine the dependence of c° on sea 
state and wind speed (Guinard, et al., 1971; Daley, 1973). The determination 
that a 0 was in some way proportional to wind speed led to the proposal for a 
satellite-borne radar, by then known as a scatterometer, to obtain oceanic 
wind and wave information (Moore and Pierson, 1966). However, at that 
time, a successful correlation of o c with wind speed was yet to be established 
and the technology with which to build a stable, accurate (±2dB absolute) 
radar was not yet available (Jones et al., 1982). 

During the 1970's, the National Aeronautic and Space 
Administration (NASA) sponsored scatterometer research and development 
which led to a technically improved aircraft scatterometer (Advance 
Applications Flight Experiment Radiometer-Scatterometer, AAFE- 
RADSCAT) and a proof-of-concept spacecraft scatterometer (Skylab, S-193). 
With the RADSCAT instrument, Jones et ah, (1977) demonstrated the 
azimuthal dependence of c°; i.e., the cosine-like dependence of o°with wind 
direction relative to the antenna. Using the radar data and corresponding 
in-situ data, empirical algorithms ("model functions") were derived to relate 
the radar backscatter to the ocean surface wind vector. In June 1978, as part of 
a suite of microwave instruments onboard Seasat, The Seasat-A Satellite 
Scatterometer (SASS) was launched. 

Three months into its planned three year mission Seasat suffered a 
failure in its power subsystem and communications with the satellite ceased; 
in November 1978 the mission was officially terminated (Pounder; 1980). Ten 
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years later the scatterometer data from this abbreviated mission are still being 
analyzed to confirm the instrument as a breakthrough in maritime 
meteorology and oceanography, but also to reveal the shortcomings in the 
scatterometer model function (Woiceshyn, et al., 1986). 

Follow-on missions to the "proof-of-concept" Seasat were already 
planned at the time of its failure. NASA had proposed a Seasat-B mission in 
the summer of 1978. Early in 1979 there was an interagency (NASA, Dept of 
Defense, Dept of Commerce) proposal for a National Oceanic Satellite System 
(NOSS, Joint Effort by NASA, NOAA, Navy, 1979 internal report). Both of 
these systems were to be "high-heritage" and "limited-operational;" i.e., they 
would carry microwave instruments nearly identical to those on Seasat and 
no new, high risk (developmental) instruments would be flown. There were 
some changes proposed for NOSS: the scatterometer was re-designed to have 
three instead of two antennas, no synthetic aperture radar would be flown 
and a large-aperture multichannel microwave radiometer was to be 
developed. NOSS grew beyond its original charter but still survived internal 
review processes until 1981 when it was cancelled. 

Reiterating the need for an improved, all-weather oceanic data base, 
the Navy proposed N-ROSS (Navy Remote Ocean Sensing System), a 
supposedly "scaled-down" version of NOSS, in April, 1981 (Honhart, 1984). 
This satellite was to be the centerpiece of a constellation of ocean observing 
satellites to be launched in the early 1990's by the US, Europe and Japan. 

t 

Funding and responsibility for the program was primarily shared by Navy 
and NASA; NOAA and the Air Force were to provide data acquisition and 
distribution support. Citing an unacceptable cost growth, the Navy cancelled 
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the program in December, 1986 (Graham, 1987; Matthews; 1987). The 
replacement launch vehicle for the N-ROSS NASA-Scatterometer (NSCAT) 
is currently planned to be the Japanese Advanced Earth Observing Satellite 
(ADEOS), tentatively scheduled for launch in 1995. Therefore, work still 
continues on the data processing algorithms including the model function. 

The European Space Agency will launch a scatterometer on its Earth 
Resources Satellite-1 in 1990. Looking to the next century, NASA scientists are 
proposing an alternate design for the Earth Observing Satellite system 
scatterometers which will use scanning instead of fixed antennas. A scanning 
system is more difficult to stabilize on orbit, but these antennas are more cost- 
effective and energy-efficient. Furthermore, they are able to operate as 
radiometers as well as radars so that it is possible to obtain corresponding 
passive microwave measurements such as atmospheric moisture parameters. 
With these changes the scatterometer will become a far more attractive 
instrument in a programmatic sense (Brown and McCandless, 1988), thereby 
increasing the chance that the US w T ill have a scatterometer in space 
sometime during the next decade. 

C ELECTROMAGNETIC SCATTERING THEORY 

While scatterometer programs have been proposed and cancelled in the 
recent past, advances in knowledge of the scattering physics have been steady 
(Rice, 1951; Barrick, 1968; Wright, 1966, 1968; Valenzuela, 1968, 1978; Fung and 
Chan, 1969; Keller and Wright, 1975; Brown, 1978; Bahar, 1981; Durden and 
Vesecky, 1985; Plant, 1986; Holliday, 1986; Donelan and Pierson, 1987). To 
date, the three most common approaches to scattering problems are physical 
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optics, small perturbation and two-scale composite surface methods 1 . Details 
of the derivations for these methods can be found in recent publications 
(Ulaby et al., 1982; Durden, 1986; Plant, 1988). Here, only the assumptions, 
results and experimental support for these theories are summarized. 

1. Fundamentals of Surface Scattering 

An expression for the NRCS (see Appendix A) is in terms of the 
electric field incident on the surface, Ej, and the field scattered from the 
surface, E s , for a given incidence angle, 0, azimuth angle, x, polarization, p, 
illumination area, A, and range, R: 



i.e., the computation of the cross-section requires specification of the scattered 
field. Fundamentally, this means that the wave equation for E s is solved 
using Green’s second vector theorem (See Appendix A, Beckmann and 
Spizzichino, 1987; Ulaby et al., 1982). This theorem states that the field at any 
point within a source-free domain bounded by a closed surface S can be 
expressed in terms of the tangential fields on the surface. 

For a surface rough in two-dimensions, the "Stratton-Chu" 
equations result as follows: 



] An alternative method is the generalized "full-wave” approach 
employed by Bahar (1981) and Brown (1978). Due to its generality, the 
approach and results are considerably more complicated and less easily 
interpreted in terms of the ocean wave spectrum. Therefore, this method has 
not been a widely adopted solution method. 




(1) 




(2a) 
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(2b) 




where a time factor of e i(0t is understood; E, H are the total electric and 
magnetic fields on the surface, K is the scalar Green's function; R is the range 
to the surface; r s is the unit vector in the scattered direction; r| 0 is the intrinsic 
impedence of the scattering medium, and k em is the electromagnetc (em) 
wavenumber (Ulaby, et al., 1982). As required by Green’s theorem, the 
scattered field is given in terms of the tangential fields, n a E and n a H. 

These equations simplify to the "Helmholtz Integral" for a surface 
rough in one dimension from which it is seen that the scattered field at 
observation point r above the surface requires the field value and its normal 
derivative at the scattering surface 



In equation (3) G is the two-dimensional Green’s function and the magnetic 
field H replaces E in the case of vertical polarized incident radiation. 

The essence of the scattering problem is the specification of the total 
field on the scattering surface. If the surface is irregular, then no general 
solution for the fields exist (Ulaby, et al., 1982) and three solution methods are 
possible: 1) the field is expanded as a series and the exact boundary conditions 
are applied (Brown, 1978; Bahar, 1981); 2) the integral equations are solved 
numerically (Durden, 1986); or 3) approximations of the boundary conditions 



,nAE are made, allowing closed form solutions to be derived as 




(3) 
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outlined below. The simplifying assumptions generally include one or more 

of the following (Beckmann and Spizichinno,1987): 

1) (often) the surface is perfectly conducting, 

2) the dimensions of the scattering elements are much larger or smaller 
than the electromagnetic wavelength, 

3) the radius of curvature of the scattering elements is much larger than 
the electromagnetic wavelength, 

4) shadowing effects from larger elements are neglected, 

5) only the far-field values are computed, 

6) multiple scattering is neglected, 

7) the density of the irregularities (per unit area or volume) is ignored, and 

8) treatment is restricted to a certain model of surface roughness. 

These assumptions al 1 relate to the reflection properties of the 
scattering surface which depend on its electrical properties and roughness. 
The roughness of a random surface is radar wavelength dependent and is 
characterized by the rms surface elevation, h, and the surface correlation 
length, /. The appropriate approximations to use for a particular problem are 
based on certain limiting values of these two parameters. 

To be "smooth," the rms height must satisfy the Fraunhofer 
condition: 



h < 



A 

32sin 6 



(4) 



where A is the electromagnetic wavelength and 0 is the incidence angle. This 
condition guarantees that there will be no more than a n/8 difference in 
phase between two reflected waves. For wavenumber k - 2n / X, the limit for 
a smooth surface is k/t < 0.3 (Ulaby et al., 1982) and the correlation length, /, is 
infinite since every point on the surface is correlated with every other point. 
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Reflections from a smooth surface are "specular” or mirror-like in that the 
angle of reflection equals the angle of incidence and very little of the radiation 
is scattered (Figure 3). Therefore, a monostatic (single antenna) radar must 
transmit at normal incidence to receive backscattered power from a smooth 
surface. 




Figure 3. Examples of Surface Scattering Patterns 
2. The Physical Optics Approximation 

The near normal condition can be met for the sea surface at 
incidence angles less than about 20° since it is possible to have wave slopes of 
similar degree to maintain the required normal incidence. In the physical 
optics, or Kirchoff approximation, the roughened sea surface is approximated 
by small planar "facets," each locally tangent to the sea surface as shown in 
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Figure 4a. Each facet specularly reflects the incident radiation so that the total 
surface field is the sum of the incident and reflected fields 

E = (1 + R)E, (5a) 

and 

^ = (1 -R)Ej-n (5b) 



where R is the Fresnel reflection coefficient of the tangent plane at a 
particular point and r is the unit vector in the scattered direction. These are 
the boundary condition approximations for the physical optics method. 
These approximations are valid when the radius of curvature of the surface 
elements is much greater than the electromagnetic wavelength and the 
surface roughness is comparable to or larger than the electromagnetic 
wavelength (see Figure 4a). 

The cross-section is then computed as the sum of the reflected energy 
from all the facets oriented normal to the radar in the illuminated area using 
the probability distribution of surface slopes. Assuming an isotropic rough 
surface and Gaussian distributed surface slopes, Barrick (1968) determined the 
cross-section to be 




sec 4 6c 



(tail 3 e/s 3 j 



( 6 ) 



where R(0) is the Fresnel reflection coefficient at normal incidence and s 2 is 
the total slope variance of the wave slopes in the x and y directions. If the 
surface is non-isotropic, then the cross-section in the upwind direction is 
written 
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(7) 



J QS 



|*(q)T 

2s u s c 



(sec 4 0) 



- tan' 1 6H l 



where s u and s c correspond to the slopes in the upwind and crosswind 
directions (Valenzuela, 1978). The slopes of importance here correspond to 
ocean wavelengths much longer than the microwavelength (to satisfy the 
criterion of small radius of curvature). The cross-section in (7) is termed 
quasispecular since it arises from only those parts of the sea surface which are 
normal to the incident wave vector. 

3. The Small Perturbation Approximation 

For a slightly rough surface and beyond about 30° incidence angle, 
comparable wave slopes are not physically realizable so that specular 
scattering is no longer valid, although scatter is still observed. At these 
greater angles it has been shown (Wright, 1966) that the response is from a 
resonance phenomenon between the em waves and the surface height 
displacements of the surface waves. This phenomenon, known as 
Bragg-scatter or Bragg-resonance, occurs when the em and ocean wavelengths 
(X, L) satisfy the equality: 



L = 



nX 

2sin0 



n = 1,2,... 



(8) 



As shown in Figure 4b, when the excess distance from the source to each 
successive crest is nX/2, the round trip phase difference is 360° and the 
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Figure 4. Schematic Drawings of the Sea Surface for a) the Kirchoff (Tangent 
Plane or "Facet") Approximation, b) the Bragg-Resonance Model, and c) the 
Two-Scale Composite Surface Model Where X cnv kenv Electromagnetic 
Wavelength and Wave number; h. Height above Mean Level; 1, Correlation 
Length; L, Xb: Ocean and Bragg Wavelengths; Longer Wave of Two-Scale 

Surface 
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reflected signals add in phase. For example, for a radar wavelength of 2.18cm 
at 0 = 40°, the Bragg wavelength is 1.7cm corresponding to an ocean capillary- 
gravity wave. 

For such small waves it would seem that the resonant signal would 
be swamped by the return from the much larger waves. Average power 
return from the larger, randomly spaced waves is proportional to the number 
of the scatterers in an illuminated area. However, in the case of resonance, 
the average power is proportional to the square of the number of scatterers. 
Therefore, the resonant effect produces a significant power return, as well as 
providing a wavelength selective tool for examining the sea surface (Ulaby, et 
al„ 1982). 

Assuming that the backscatter is from very short, small amplitude 
waves, the cross-section is derived using a perturbation approach in which 
the mean scattering surface is flat and small wave perturbations are 
superimposed (Figure 4b). "Small" perturbations are defined to have rms 
surface slopes which are sufficiently gentle (V 2 h/l < 0.3) and surface 
displacements small compared to the em wavelength (k h < 0.3). 2 The total 
field on the surface is written as the sum of the incident and reflected fields 
(specular reflection from the flat mean surface) plus a small scattered 
component from the small amplitude roughness. 



2 Chen and Fung (1988) suggest "that a more appropriate description for 
the range of validity of this model is a two-dimensional space with kh and kl 
as the axes"; i.e., in addition to small height and slope the horizontal 
correlation length must also be large compared to the electromagnetic 
wavenumber. 
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The rough surface is considered random within a particular L by L 
area but this random surface repeats itself in distance, L; i.e., is periodic. 
Flouquet's theorem is then invoked which states that a scattered field from a 
periodic surface is itself periodic. Therefore, it is appropriate to write the 
scattered field as a harmonic series expansion: 



wavenumber component in the z direction. That is, the scattered field is 
periodic in x and y and propagates upward. 

The relationship between the coefficients, B mn , for the different field 
components is found from Maxwell's divergence equation V- E = 0; .i.e., there 
are no sources in the region. The boundary condition of perfect conductivity 
(the tangential fields on the surface are zero) provides the additional two 
equations needed to solve for the three coefficients. Since it was assumed that 
k h was small all exponentials involving z are expanded in Taylor's series of 
k h about h = 0 and the coefficients are expanded as increasing powers of the 
surface elevation. For first-order Bragg scattering only the first terms are 
retained. The dominant effect of the first order expansion is to give rise to a 
z-component in the scattered field. Higher order expansions (Valenzuela, 
1968) produce cross-polarization effects. 

The radar cross-section for resonant scatter from a dielectric 
(nonconducting) surface is derived to be (Wright, 1968) 



E s = 2X „E{P + m ' n ' z ) 



(9) 



where E(p + m,n;z) = e 




e L e 



, p = 2k / X. sinG and k z is the 



o°{9) = 16 nk 4 cos 4 0g t} 2x ¥(2ks\nd,0) 



( 10 ) 
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where 



Shh ~ , 

cos 6 + (e, - sin 6 )? 



(Ha) 




(e r - l)je,(l + sin 2 0) - sin 2 



(lib) 



and ^(kx, ky) is the two-dimensional ocean wavenumber spectrum. The 
complex relative permittivity of the ocean, Er, is given by e r = e/e 0 - io/£ 0 (0, in 
which the real part represents the capacity of the medium to store electric 
energy (c.f. a vacuum for which e G = 8.85X10‘ 12 Farad/m). The imaginary part 
represents losses by the medium through the conductivity, o, at frequency, co. 
The subscripts refer to the transmit/receive polarizations of the signal: h- 
pol/h-pol and v-pol/v-pol. For a perfect conductors = ■», e r = -i «*, ghh= 1 and 
gw = (1 + sin 2 0) / cos 4 0. Therefore, resonant scattering of h-pol radiation 
from a good conductor is more sensitive to incidence angle than is v-pol. 

In his classic wavetank experiment, Wright (1966) demonstrated that 
some of the observed characteristics of backscatter can be explained by this 
"first-order" Bragg theory. As shown in Figure 5, the measured backscattering 
and the polarization ratio (o V v/Chh) as functions of incidence angle agree 
rather well with the theoretically computed curves. Wright hypothesized 
that the differences from the theoretical cross-sections could be due to 
oversimplifying assumptions regarding the electric fields near the surface 
when no waves were present. 



25 



•10 



Vtftical 





Backscaitenng cross section as a function of incidence angle for capillary waves on 
tap water. Left; the scattering section calculated from theory compared with observations. Right' 
the ratio of vertically to horizontally polarized backscatter, also compared with theory (from 
Wnght. 1966). 



Figure 5. Evidence of First-Order Bragg Scattering 
(Stewart, 1985 after Wright, 1966) 

Further evidence for Bragg-scattering is found via analysis of the 
Doppler spectrum of the scattered signal. The Doppler spectrum is plotted as 
the spectral density of radar returns against Doppler frequency shifts "induced 
in backscattered microwaves due to a surface moving with a line-of-sight 
velocity toward or away from the radar" (Plant and Keller, 1989). If there are 

t 

no background currents and Bragg scattering is the appropriate mechanism, 
then the Doppler spectra will contain peaks at those frequency shifts 
corresponding to plus and minus the Bragg wave frequency. 
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Plant and Keller (1989) provide evidence of the Bragg scattering 
mechanism in field data. Doppler spectra computed from an L-band (16cm) 
scatterometer positioned at the end of the Coastal Engineering Research 
Center pier near Duck, North Carolina are shown in Figure 6. As predicted, 
spectral peaks appear at the Bragg wave frequency of 3.35Hz for the 
microwave parameters used in this experiment (Figure 6a). The wind 
direction for this spectrum is 210°, i.e., the offshore direction, and the wind 
speed is only 4.3m/s; hence, the dominant ocean wavelength will be 
correspondingly short. On the other hand, the peaks disappear for a spectrum 
corresponding to the onshore direction (Figure 6b). The longer waves have 
orbital velocities which contribute to the Doppler shift, thereby broadening 
the Doppler spectrum. It appears then that the first-order theory is valid until 
the dominant wind wavelength slightly exceeds the Bragg wavelength (Plant 
and Wright, 1977). At this point, the k em h<t<tl and the mean surface is no 
longer flat. 

4. The Two-Scale Composite Surface Approximation 

These limitations of the pure first-order theory led Wright (1968) and 
Bass et al., (1968) to propose independently a composite surface (two-scale) 
theory in which the surface is envisioned to consist of a large scale part which 
has small curvature and a small-scale part which backscatters according to 
first-order theory. The sea surface is broken up into an infinite number of 
slightly rough patches and the net backscattered power is computed as the 
summation of the power from a single patch distributed over 0 the slopes of 
the dominant wind waves of the ocean (Figure 4c). The wavenumber 



27 



spectrum is decomposed into two regions corresponding to the longer tilting 
waves and the shorter Bragg-scattering waves. 



Field Data Two-Scale Model 





(a) Offshore Wind Direction 





(b) Onshore Wind Direction 



Figure 6. Evidence for First-Order Bragg Scattering and Validity of Two-Scale 
Model (a) (b) Left, Doppler Spectra from Field Data; (a) (b) Right, Doppler 

Spectra from Two-Scale model 
(adapted from Plant and Keller, 1989) 

The radar cross-section is the sum of two terms: the quasispecular 

contribution from the tilted surfaces and the Bragg contribution from the 

small amplitude waves. 

o°(0) = e~ il<>, 1 > CT’ QS + | a i (d,y,8)p(tan y,tan<5)dtan >dtan<5 (12) 

where 
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<7b(0/7,<5) = 16^ 4 cos 4 0|g(0,y,<5)|y/ s (2):sin(0 + y),2kcos(6+ y)sin<5) (13) 



is the Bragg-scatter cross-section for a slightly rough surface tilted by y in the 
x-direction and 8 in the y-direction. 

In (13), Ys is the short wavelength portion of the wave spectrum. 
The quasispecular term in (12) is now reduced by a factor exp(-4k 2 </i 2 >) due 
to the small waves which make the specular "facets" less smooth, and the 
Bragg contribution is altered by the presence of the longer waves which tilt 
the Bragg waves in and out of the plane of incidence (Valenzuela, 1968). A 
separation wavenumber must be defined between the two assumed scales for 
the composite surface to be applied correctly. If such a "spectral gap" exists, 
then the scales are effectively decorrelated, allowing the linear summation 
of the two types of backscattered power. 

The relationships between backscatter and incidence angle 
corresponding to the three types of scattering models discussed are illustrated 
in Figure 7. Composite theory compares well with observed values of a 0 for 
vertical polarization, with an assumed ocean wave spectral form proportional 
to k' 4 , but the curves diverge from the observed values for h-pol at higher 
incidence angles. Assuming higher, somewhat unrealistic, background slopes 
(10° and 20°) helps to fit the data but the remaining differences imply that 
other scattering mechanisms may become important for h-pol scattering at 
grazing angles (Valenzuela, 1978; Donelan and Pierson, 1987). Better 
agreement (not shown) between observations and theory are found at lower 
microwave frequencies (.438 Ghz, L-band) than higher frequencies (4.4 Ghz, 
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Com pan son of measured arid theoretic*] cross sections of the ocean for a radar fre- 
quency of 4 5 GHz Topr vertical polarization. Bottom: honzontal polanzation The vanation in 
the observations is due pn manly to vanation in wind speed at the sea surface (from Valenzuela, 
197*), 



Figure 7. Backscatter Versus Incidence Angle for Three Scattering Regimes 
(from Stewart, 1985 after Valenzuela, 1978) 
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C-band) (Daley, 1973). Finally, composite surface theory explains the widened 
spectral peaks observed in the Doppler spectra as being due to the advection of 
the small-scale structure by the orbital velocity of the longer waves (Figure 
6b) (Keller and Wright, 1977). 

D. SPECIFICATION OF THE ROUGH SEA SURFACE 

The expressions for the radar cross-section in the previous section show 
that the scattering cross-section is linearly dependent on the two-dimensional 
ocean wavenumber spectrum. To examine pure Bragg scattering with the 
first order theory, the spectral density for a small range of w T avenumbers near 
the Bragg wavenumber is required (the high frequency tail of the full 
spectrum). Two scale composite surface theory also requires the joint 
probability density function of the upwind and crosswdnd slopes or, 
alternatively, the low r er frequency gravity wave part of the wavenumber 
spectrum. Although the correspondence is good between the observed and 
computed backscatter, the absolute accuracies of the Bragg and two-scale 
formulations cannot be unambiguously established since a universally 
acccepted functional form of the tw T o-dimensional w'avenumber spectrum has 
not yet been derived. Thus, present model function research is necessarily 
semi-empirical and semi-theoretical: wave and radar data are used in a 
feedback loop with w T ave theory to derive new forms for the wavenumber 
spectrum and to determine the physical variables of importance to the 
development of the ocean weaves, thence backscatter. It is at this point that 
the transition is made from a pure electromagnetic scattering problem to an 
ocean wave dynamics problem. 
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1. Ocean Wave Spectra 

The ocean surface is an evolving, random rough surface composed 
of waves which are generally categorized by their generating or restorative 
forces. The surface waves visible to the observer are most often the result of 
the driving force of the wind and the restorative forces of surface tension 
(capillary waves) or gravity (wind generated gravity waves). These waves 
span wavelengths of millimeters to hundreds of meters (wavenumbers of 
1000 to .001 m' 1 ) and have periods of less than 0(1 0s) (intrinsic frequencies of 
greater than 0.1 s' 1 ). Wave triads or quartets (capillary, gravity waves, 
respectively) can interact nonlinearly, with the result of an energy transfer 
among the waves, provided the wavenumbers of the individual waves satisfy 
certain resonance conditions (Hasselmann, 1962). Energy is removed from 
the wave system by wavebreaking, the generation of "parasitic" capillaries on 
the steeper gravity waves, or frictional dissipation at the highest frequencies 
(Hasselmann, et al., 1973). 

All of these competing generation and dissipation mechanisms lead 
to a random sea surface which is best described by its overall statistics, as 
formalized by spectral (Fourier) methods; i.e., as a linear superposition of an 
infinite number of sinusoidal components (Pierson, et al., p.24, 1955). This 
requirement is satisfied if the components are of small amplitude and slope 
and can be approximated by a sinusoidal profile. Under these circumstances 
the Fourier transform of the surface elevation covariance is computed as a 

t 

function of wavenumber and frequency, with position and time as 
parameters. This "power spectrum" is the spectral "signature" of the sea 
surface. Donelan, et al., (1985) point out that the higher-order effects of 
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deviation of the wave profile from a pure sinusoid, nonlinear interactions 
and viscous dissipation will prohibit the signature from being permanent, but 
to first order it is an accurate representation of the wave field. 

Assuming the wave field is spatially homogeneous and statistically 
stationary, then the ergodic theory applies: space and time averages from a 
single time series (realization) are equal to ensemble averages across all 
possible realizations (Bendat and Piersol, 1986, p. 144). This random ocean 
surface is then characterized by the covariance, p(r, t), of the surface 
displacement, £, at points separated by a distance r(x,y) and time, t 

p(r,f)=<s(x,f 0 )f(x + r,f 0 + f)> (14) 

where the angle brackets denote an ensemble average. The covariance 
quantifies the decorrelation scales of the sea surface in space and time. There 
are then three spectra of varying dimension that can be derived from this 
covariance (Kinsman, 1965). 

The three dimensional (k x , k y< to) wavenumber-frequency spectrum 
is defined as the Fourier transform of the covariance (Wiener-Kinchin 
theorem, LeBlond and Mysak, 1978, p. 311) 

S(K, oj) = j P( r ' ty^'^drdt (15) 

and inversely, 

p(r, t ) = J J Js(k, coy^'-^dVda (16) 

At r=0 and t=0 equation (15) gives the variance < £ 2 > of the surface elevation, 
or mean-square elevation. The calculation of S(k,co) requires simultaneous 
measurements of surface displacement in time and space, an infrequently 
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performed task which requires an array of wave staffs and intensive data 
processing (Donelan, et al., 1985). 

The two-dimensional wavenumber spectrum required for the NRCS 
calculation is the integral of the three-dimensional spectrum over all 
frequencies 

T'(k) = js(k,co)dco (17) 

At t=0, the covariance is 

p{ r) = p(r,0) = Jj (18) 



then 

4,(k)= (^f!! p(r,e " k,>rfr <19) 

This spectrum is also difficult to measure since it requires surface elevation as 
a function of space, an inherently more difficult in-situ measurement to 
make which involves the use of either stereophotography (Shemdin, et al., 
1988), laser optic techniques (Tang and Shemdin, 1983; Keller and Gotwals, 
1983) or a spatial array of sensors (Donelan, et al., 1985). 

It is much more common to have only a single wave staff acquiring 
data in time. In this case, a one-dimensional frequency spectrum results 
which is independent of wavenumber magnitude or direction 

Q(co) = fjs{k,co)dk ( 20 ) 
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p(t) = p(0,t) = JJ Cl{co)e ib *da> 



( 21 ) 



and 



n(a) = ^ J p{,)e ‘“ dl 



( 22 ) 



or 




(23) 



since Q(cd) is symmetric about co=0. 

An approximate wavenumber spectrum can be determined from the 
frequency spectrum using the dispersion relation, co=co(k). The wavenumber 
spectrum is often expressed in polar coordinates as the product of a 
wavenumber magnitude and direction spectrum 



where k = (k cos (p ,k sin (p),k = |k|, and cp is the wavenumber bearing (Pierson, 
1955). The wavenumber spectrum is obtained from the frequency spectrum 
using the dispersion relation from which 



(LeBlond and Mysak, 1978). This procedure cannot be used to determine the 
large wavenumber portion of the spectrum since the intrinsic frequency of 
the shortest (capillary) waves, co 0 , is Doppler shifted by the orbital velocity of 
the longer waves on which they ride. The orbital velocity is essentially a 



v F(k) = v P(]t,<p) = SOt)0(<p) 



(24) 




(25) 
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slowly varying current of velocity U(x,t) so that the frequency of the short 
waves is given by 

a>(k) = ty 0 (^) + k-U(x / f) (26) 

That is, the relationship between frequency and wavenumber is no longer 
unique (Atakturk and Katsaros, 1987). 

Higher order statistics can be obtained from the wavenumber 
spectrum. For example, the mean-square slope of the waves is given by the 
second moment of the wavenumber spectrum (Cox and Munk, 1954) 

s 7 =jk 2 '¥{k)dk (27) 

The slope spectrum is often used in the study of the high frequency portion of 
the wave spectrum since higher wave numbers are weighted more heavily. 
Cox and Munk (1954) found that the mean-square slope is wind speed 
dependent. For a clean surface and windspeed at 12.5m in m/s, 

s 2 =1 O' 3 (3 + 5.1 2U) (28) 

and 

s 2 = 1CT 3 (8 + 1.5617) (29) 

for a slick-covered surface. In the Cox and Munk study, the slick was 
artificially generated and composed of diesel fuel and fish oil; the observed 
slopes may differ with typical concentrations of natural films. 
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In polar coordinates, with the x-axis pointed upwind, the 
mean-square slopes in the upwind /downwind and crosswind directions are 



s x = JJ ( kcos(pf x i'(k / (p)dkd(p (30) 

— O© 

and 

Sy=jJ ( k sin (pf'¥(k,(p)dkd(p (31) 



These quantities are of particular interest in the present study since they are 
required for the computation of the quasispecular and Bragg scattering terms 
of the two-scale model and for the expressions describing the hydrodynamic 
modulation of the Bragg waves by the longer waves. 

2. The Spectral Transport Equation 

The adoption of spectral forms in wave analyses is a statistical means 
of revealing basic structures in the random wave field. Underlying these 
statistics are, of course, the physical mechanisms responsible for the waves. 
The individual wave components generate, propagate and decay at different 
rates so that the spectral shape changes as the wave field develops. LeBlond 
and Mysak (1978, p. 322) then pose the question 

...as to why spectral (and other statistical properties) take the shapes they 
do at various stages in the life of a wave field. Is it possible to derive 
predictive (as opposed to purely empirical or phenomemological) 
theories for the statistical properties of varying wave fields? More 
specifically, given the statistics of the forcing functions (such as wind, the 
atmospheric pressure), the boundary and initial conditions and the 
coefficients of the governing differential equations, each of which will, 
in general, have a random component, can we compute the solution and 
its statistics in terms of the statistics of the forcing? 
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The "radiative transport theory" of the evolution of spectra, which 
arises in the area of stochastic fluid dynamics (Monin and Yaglom, 1965,1967), 
provides the mathematical formalism to answer these questions. The 
application of the theory (Hasselmann, et al., 1973) rests on the assumption 
that the time and space scales of spectral evolution are much longer than the 
longest resolvable period and wavelength in a data record. This assumption 
allows the evolution process to be weakly nonstationary. It is found in the 
absence of currents that the ocean wavenumber spectral density is conserved 
along a wavenumber path except for certain exchange processes, denoted by L. 



These exchange processes include the energy input by the wind Qj; 
Qnb the redistribution of energy due to nonlinear interactions and the 
generation of parasitic capillaries on the steeper gravity waves; and Qd, the 
dissipation of energy by wave breaking and frictional dissipation,. The group 
velocity of the waves, Cg = dco/ dk, is the velocity at which the wave energy 
propagates. 

A more general expression which allows the waves to be 
superimposed on a space- and time-dependent mean flow, U(^,t) is in terms 
of the wave action density,N(k;x,t) 




(32) 



or 




(33) 



dN 

dt 



(34) 
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or, 



3n r 
at + ' 




(35) 



where N = 'F(k) / 0) o . The frequency is now defined by equation (26) where 
O) 0 is the intrinsic frequency of the wave in a coordinate system which is 
moving with the current, U. In the presence of such a current, the wave 
action is conserved, whereas the spectral density is not. 

A visual summary of these equations is given in Figure 8 which 
shows a typical frequency spectrum and the associated source terms 
(Hasselmann, et al., 1973). Wind energy is primarily input at and above the 
frequency of the spectral peak and removed by dissipation at frequencies 
greater than approximately three times the peak frequency. The nonlinear 
terms transfer energy from the spectral peak to both higher and lower 
frequencies, but at the higher frequencies it is cancelled by the dissipative 
effects. Different stages of the wave spectrum evolution can be evaluated by 
examining different balances among the terms in the radiative transfer 
equations (32) and (34). When d'F/dt >0, the wave spectrum is in a growth 
stage, a fully arisen sea is represented by dTVdt=0 and the spectrum is in a 
decaying state when d4Vdt<0. 

3. The Fully Arisen Sea, dTVdt = 0 

The wave spectrum cannot grow without limit. Instead, for a given 
constant wind speed, it will reach some maximum state dictated by the length 
of time the wind blows (duration) and the distance over wh'ifch the wind 
blows (fetch). At this point d4 / /dt=0, which implies that the source and 
dissipation terms are in equilibrium at all frequencies 
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Figure 8. The Mean JONSWAP Spectrum and Source Functions L'(co) 

/ 

together with the Computed Nonlinear Energy Transfer 
(from LeBlond and Mysak, 1978, after Hasselmann, et al., 1973) 

Q, +Qni + Qd=o ( 36 ) 

At wavenumbers below the wavenumber of the spectral peak, k < k p , the 
phase speed of the long waves is greater than the wind speed (for a fully 
developed sea) so that the waves lose some of their energy to the atmosphere 
only to have it restored by nonlinear transfer of energy from the spectral peak. 
These long waves are still in a growth stage and can accept more energy 
without breaking. Near the spectral peak, k = k p , wind input is transferred to 
higher and lower frequencies by nonlinear interactions. Above the spectral 
peak but below the wavenumber where viscous effects are predominant, 
1.5k p < k < 3k p , (Leykin and Rozenberg, 1984; Donelan, et al., 1985), input 
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from the wind and nonlinear interactions are balanced by dissipation 
processes, including wave breaking and the formation of "parasitic" 
capillaries ahead of the steep wave crests of the longer primary waves. 
Phillips (1958) identified this part of the equilibrium spectrum as the 
"saturation" range: the spectral density of these wavenumbers is limited by 
the onset of dissipation processes. 

Observed spectra, appropriately normalized, exhibit remarkable 
similarity with the following characteristics. The wavenumber of the spectral 
peak decreases with increasing fetch and increasing wind speed and the 
spectral density increases with increasing fetch and wind speed (Figure 9a) 




Figure 9a. The Evolution of Surface Wave Spectra with Fetch, from the 
JONSWAP Observations. The Fetch Increases with the Number labelling the 

Spectral Peak. 

(From LeBlond and Mysak, 1978, after Hasselmann, et al., 1973) 
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(Kitaigorodskii, 1962; Hasselmann, et al., 1973). There is very little energy 
below the well-defined spectral peak, while above the peak, the spectral 
density decreases more slowly (Figure 9b) (Donelan, et al., 1985). Slope spectra 
(Cox, 1958; Tang and Shemdin, 1983; Shemdin and Hwang, 1988) also display 
characteristic features. These similarities have led researchers to compute 
parametric fits to the spectral shape. 







Figure 9b. Normalized Frequency Spectra Grouped into Classes by U^Cp. The 
Vertical Bars at the Top of the Figure are an Estimate of the 907 c Confidence 
Limits Based on the Standard Error of the Mean. 

(from Donelan, et al., 1985) 
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The analytic fitting functions of the earlier models (before 1977) are 
based on Phillips' concept of a saturation range defined above. Phillips (1958) 
argued that wavebreaking in the saturation range was a dynamic instability in 
which the local acceleration of the fluid parcel exceeded the gravitational 
acceleration, g. Therefore, the spectral density was defined as a function only 
of frequency (or wavenumber) and gravity. Since the wavenumber spectral 
density has units of L 4 , by dimensional considerations we have 

S{k) = pk~* (37) 



or, in terms of frequency 



Q(co) = ag 7 0) 5 



(38) 



where a and P are universal constants. The spectral slope is either k^or co~ 5 
and independent of wind speed in the saturation range. Equation (38) is used 
as the "kernel" of the spectral model proposed by Pierson and Moskowitz 
(1964) and the Joint North Sea Waves Project (JONSWAP) model proposed by 
Hasselmann, et ah, (1973). 

Phillips (1977) modified the form of the saturation range to include 
the possibility of "microscale" breaking of the short waves (O(lOcm)) in the 
presence of a significant wind drift layer, produced by the tangential 
component of the wind stress (Banner and Phillips, 1974). Banner and 
Phillips proposed that the drift current augments the phase speed so that 
waves whose phase speeds are a small multiple of the surface drift break 
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